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FORMULAS FOR nth ORDER DERIVATIVES OF HYPERBOLIC 
AND TRIGONOMETRIC FUNCTIONS 
by Edwin G. Wintucky 
Lewis Research Center 

SUMMARY 

Formulas are derived and presented in the form of finite series for derivatives of 
any order of the hyperbolic cotai^ent, tangent, cosecant, and secant. The coefficients 
have a simple recursive property which facilitates their computation. A representative 
derivation and proof by mathematical induction are given for the hyperbolic cotangent. 
An application of the formulas to the evaluation of certain Fourier sine and cosine inte- 
grals is demonstrated. A method for obtaining formulas for the derivatives of the cor- 
respondii^ trigonometric functions is also presented. 


INTRODUCTION 

A method for evaluatii^ cubic lattice sums which arise in the theory of magnetism 
has recently been developed (ref. 1). This method used infinite series expansions which 
are partially summed by means of the Laplace transform (ref. 2) and result in certain 
Fourier sine and cosine integrals. 

These integrals and several other Fourier sine and cosine integrals are presented 
in standard tables of integrals (refs. 3 to 5 and references therein) in terms of deriva- 
tives of hyperbolic functions such as ctnh and tanh. The general case involves the 
derivative of nth order. Examples of these are (ref. 1) 



dx= (-l)*^-^ 


2n 


da 


2n 


— ctnh na. 
2 



( 1 ) 


(a > 0) 



ctnh n-a - — \ (2) 

2a j 

(a > 0) 

In reference 6, the transforms of these integrals are listed as functions of Riemann 
zeta functions in the form of infinite series, which are inconvenient to evaluate. Formu- 
las for the higher derivatives in equations (1) and (2) do not appear in any of the standard 
references on mathematical tables and formulas (refs. 7 to 11) or treatises on applied 
mathematics (refs. 12 to 18). Direct computation of a higher order derivative becomes 
inconvenient in the absence of a general formula. Furthermore, in the problem men- 
tioned previously, the integrals in equations (1) and (2) appear as the nth terms in 
infinite series. These series are more easily handled with the nth term expressed in a 
more analytical form. 

In this report, general formulas are derived which give the derivatives of the 
hyperbolic cotangent to any order in the form of finite series. Numbers are defined for 
the coefficients of the series which have a simple recursive property and are easily 
calculated; thus, the formulas are i«.rticularly suitable for numerical evaluation by 
desk -top computers with built-in programs for hyperbolic and trigonometric functions. 
Parallel formulas are also presented for the hyperbolic functions tanh, sech, and csch 
and for the trigonometric functions ctn, tan, sec, and esc. A representative induction 
proof for the formulas is given in the appendix. 

DERIVATION OF FORMULAS 

The formula for the derivative of arbitrary order of the hyperbolic cotangent (ctnh) 
is derived as follows. 

Let = ctnh y = u'u”^ where u = sinh y and u’ = cosh y. Then u” = u and 
2 ^2 

(u’) = 1 + u . For successively higher derivatives, where = (d/dy)Aj^ carefully 
rearranging terms in the following way makes it possible to discern a recursive pattern 
and thus write down the general term, which will subsequently be proved: 



A^ = -u"^ 
Ag = (2:)u»u“^ 


x^n+1 cos ax ^ 

X 1 j 2 

e - 1 da 
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Ag = -(3!)u"^ - (2 ^)u“2 
A 4 = (4!)u’u"^ + (2^)(2)!u'u"^ 

Ag = -(5l)u'® - (2^ + 4^)(3l)u'^ - (22)2 u"^ 

Ag = (6:)u’u“'^ + (2^ + 4^)(4:)u’u"^ + (22)2(2:)u’u"^ 

Arj = -(7!)u“® - (2^ + 4^ + 6^)(5!)u“® - [(22)^ + 42(£2 + 42^^ (3.)u"4 _ (22)3u' 


j 

Ag = (8!)u’u"® + ^ (6!)u’u 


(2^3)^ (2^i)^ (4:)u'u"^ + (22)3{2 !)u’u 


^2=0 ij=0 


Let us define numbers such that 


W2n, 0 = 1 


^2n, 1" 

m=0 


^2n,k~ S 

m-0 


2m,k-l 


= (2n)"w2„ k-l + W2(„_i) k 

The series for k = 1 represented by j is well known and has the sum 

2n(n + l)(2n + l)/3 (ref . 19). The derivatives of Aq can be written in terms of the 
W2n,k and we have, for example, 


Ag = (8!)u’u"® + Wg ^(6!)u’u“' + 2(4!)u’u''' + W2 3(2l)u’u 




*hi»n 
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The even derivatives for arbitrary n can then be written 
A2„= (2n)!u'u-(2“+« + j [2(n - l)]!u'u-<2"-« 

+ «^2(n-2), 2 [2<“ - 2)] M'u-<2“-3) * ... „.2(4 :)u’u-5h. W 2_„.i(2!)a'u-- 

= Z '^2(n-k),k[2(" - k)]wu-2(“-«-2 (6) 


The odd derivatives can be similarly written: 


k=0 

In terms of hyperbolic functions, 

j2n 


A2n+1 = ,,[2(1 - k) + l] :u-2<"-“+l) 


2n 


dy 


j2n+l 


n-1 

ctnli y = ctnh y^^ ^2{n-k) " k)]l(csch 


(n^l) 


dy 


2n+l 


n 

ctnh y = -2] W 2 ,„.„„ 1 )^ „ [2(n - k) t l] ! (csch y)2(n-k+D 


k=0 


(n>0) 


(7) 


( 8 ) 


(9) 


Since 


d^ 1^ (-1)^! 
dy^ y y^-*-! 


( 10 ) 


The Fourier sine and cosine integrals in equations (1) and (2) can be written explicitly as 
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/ 


sin ax , _ (-1)” _2n+l d‘ 


2ti 


e^-1 


dx = fctnh y 

2 dy^n V y / 


= ;r^^+^|ctnhy^ J2(n -k)]:(cschy)2("-l^) - i^L , (H) 


c 

/ 


x2n+l ^ ^ (_i)t^+l ^2n+2 


- 1 



'»'2(n-kH.l),k[2(“ -« * 


- \ (12) 

2n+2 

y 

where y = Tra. These integrals are thus easily and conveniently evaluated for any n to 
any degree of accuracy. Equations (11) and (12) have been checked numerically for 
n = 0 to 5 and a = 1 to 5. 

Formulas for the higher derivatives of tanh, sech, and csch, which may be derived 
in a similar way, are tabulated in the next section. A method is also described for ob- 
taining the higher derivatives of the corresponding trigonometric functions from the 
formulas for the hyperbolic functions. 

The coefficients in the derivatives of sech, csch, sec, and esc consist of the sums 
of odd numbers, where 


W 


2n+l, 0 


= 1 


W, 


2n+l, 1 


(2m + 1)" 


n 


W 


2n+l,k“ ^2m+l,k-l 


= (2n4-l)V2„^lk-l + W2„-l,k 


(13) 
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A representative proof by mathematical induction for the formulas is given in the 
appendix. 

Tables for the numbers and are most conveniently generated 

k 

using W2,k=2^^ '^4,k=^22<>'+ni), ,^ = for 

2 

even numbers and Wj^ k “ ^ ^2n+l k ~ ''^2n+l k-1 ^2n-l k 

numbers. 


TABULATION OF HIGHER DERIVATIVES 

All the formulas presented in this section may be derived in the manner outlined for 
the hyperbolic cotangent in the previous section, the formula for which is repeated here 
for the sake of completeness. 


Hyperbolic Functions 


j2n 


n- 


<3y 


2n 


ctnh y = ctnh j^[2(n - k)]l(csch (14) 


j2n+l 


n 


dy 


2n+l 


ctnh y = + l] '(^sch y)2(n-k+l) (15) 


;i2n 


n-1 


dy 


2n 


tanh y = tanh y 2(n - k)]!(sech 


(16) 


q2n+l 


tanh y - k) + l] !(sech y)2(n-k+l) (17) 


dy^"+^ k=0 


csch y W2(„.y^j [2(n - k)] l(csch y)2<-«+l 

dy'*'^ k=0 


(18) 


I 



I 


j2n+l 


dy' 


2n+l 


a 

csch y = - ctnh y ^^ W2(n-k)+l,k[^^*^ - k) + l] !(csch y)2(n-k)+l (jg) 


k=0 


42n 


dy 


2n 


sech y “ k)]t(sech y)2(n-k)+l 


( 20 ) 


j2n+l 


dy 


2n+l 


n 

sech y = - tan y^ j^[2(n - k) + l] !(sech y)2(n-k)+l (gl) 


k=0 


Trigonometric Functions 


Analogous formulas for the corresponding circular functions can be simply obtained 
by making the following substitutions: 


ctnh y = ictn iy 1 


tanh y = - itan iy 
csch y = icsc iy 
sech y = sec iy 

d^ ^.m d"^ 

dy"' d(iy)"' 


(22) 


Then, for example 


n-1 


ctn y = ctn y^ k[2(n - k)] l(csc y)2(n~k) (23) 

dy2n - - 


k=0 


(n>l) 
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CONCLUDING REMARKS 


Previous untabulated formulas for the derivatives to any order of certain hyperbolic 
and tr^onometric functions have been derived and presented. These formulas are given 
in the form of finite series, the coefficients of which have a simple recursive property 
and thus are easily calculated. At least one use for these formulas is in the evaluation 
of Fourier sine and cosine integrals such as those mentioned in the report. 

As n increases, the coefficients become very large. In actual applications, how- 
ever, there may be multiplying factors which may somewhat offset this trend. Such is 
the case for the problem in the theory of magnetism referred to in the INTRODUCTION. 
For example, factors multiplying the integrals in equations (1) and (2) are of the form 

[(2n)!]^ 

[22’^(n:)2]^ 


For general reference purposes, a complete listing of the 16 formulas for the nth 
order derivatives is not necessary. The four formulas in equations (15), (17), (18), and 
(20), together with a definition of the numbers are sufficient. The other formu- 

las are then easily gotten by either a single differentiation or simple substitution using 
equation (22). 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, May 4, 1971, 

129-02. 
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APPENDIX - PROOF BY MATHEMATICAL INDUCTION OF FORMULAS 
FOR DERIVATIVES OF HYPERBOLIC COTANGENT 


A detailed proof by mathematical induction of the formula for the odd derivatives of 

ctnh y, d^^^^^cinh y/dy^’^'*'^, is given here to verify its validity for all n. The formula 
2ri 2n 

for d "^ctuh y/dy ” is consequently also verified. Similar proofs by mathematical induc- 
tion can be constructed for the tanh, sech, and csch formulas. 

The derivatives of d^*^‘*'^ctnh y/dy^*^'*'^ for n = 0, 1 are 


d . . ,2 -2 

— ctnh y = - csch y = - u 

dy 


(Al) 


ctnh y = - (3'.)csch^ y - (2^)csch^ y 


dy‘ 


(A2) 


= - (3I)u"‘^ - (2^)u"^ 


where u = sinh y. These cases are readily verified by direct calculation. 
Assume the formula true for n = m - 1. Then 


j2m-l 


m-1 


dy 


2m-l 


ctnh y . -V J2(m - k) - 


(A3) 


By direct differentiation, 


j2m 


m-1 


dy 


2 m 


ctnh y = V j^[2(m - k)] lu.u-2(-k)-l 


(A4) 
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,2m+l 


2m+l 


dy 


. ctnh y = - V " k) + l]ui + u2)u-2(“-J^+D 


L '^2(m-k),k[2(m -k)]:u' 


m-1 

= - ® 2 (m-k),k[ 2 (” - I') 

k=0 


■ £ «'2(m-k),k{[2<”> - k + l)]l - [2(m - 


-k) 


Consider the first sum. 




- «'2m,0<2"’ + I):--''"*'* * ^W2,^-k),k[2(“ ' « + l] ’.u-^fk-k-^D 


In the second sum, 


[2(m - k) + l] : - [2(m - k)] ! = [2(m - k)j^|^2(m - k) - l] ! 


The second sum is then 


2(m-k) 


(A5) 


(A6) 


(A7) 
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m-1 

2^ [2(m -k)]V 2 ,„_y ,^[2(m - k) - l]..u-2(”-« 
k=0 


m-2 
k“ 0 


[2(m - k + l)]S(m-k+l),k-l[2(”“ ‘ « + ^ Wj 


(A8) 


where the dummy index k is replaced by k - 1 and Wo = (2^)Wo i by definition. 

^ j m z j m “ X 

Recombining the two series gives 


d 


2m+l 


dy 


2m+l 


ctnh y = - W 2 j^ Q(2m + 


{[2(m - k) + 1^., J[2(m - k) + l] 


-Wo ^u' 
^ , m 


(A9) 


By the definition of the numbers W 


2n,k, 


'^2m+l,0“ ''^2m,0" ^ 


and 


[2(m - k) * l]S(m-k,l),k-l - *2(m-k),k 

= [2(m - k+ l)]S(m-k+l),k-l + 5Z(2')Si,k-l 


m-li 


m-k+1 


£ (2i)^,,k-l 


= w, 


2(m-k+l),k 


(AlO) 
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Then 


j2m+l 

-2 ctnh y = 

dy2“+l 


'"'2(m+l),0 


(2m + l)!u 


-2(m+l) 


m-1 


W, 


2(m-k+l), k 


k=0 


[2(m - k) + l] tu-2(m-k+l) ^ 


u 


-2 


m 


k=0 


W, 


2(m-k+l), k 


[2(m - k) + 


(All) 
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